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Abstract 

We check the large N transition proposed by Aganagic and Vafa and further studied 
by Diaconescu, Florea and Grassi using the large N limit of the corresponding Chern- 
Simons two matrix model. We find the spectral curve and calculate the genus zero free 
energy. 
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1 Introduction 



The conifold transition is a duality between open and closed strings. In the topological 
A-model, this is a duality between the open string model on the deformed conifold and 
the closed string model on the resolved conifold. This was originally studied by taking 
the partition function of large N Chern-Simons (CS) theory on expanded in a 't 
Hooft limit and then summing over the holes to get a closed string theory 121 • A 
worldsheet proof of this duality has since been provided [Sj and generalized in 

It is interesting to extend this duality. One possible generalization is the local 
conifold, a Calabi-Yau (CY) that locally looks like a conifold but has an extra vertex 
in its toric diagram. The large N transition for this type of geometry was conjectured 
in |B| and further investigated in jB] . The all genus partition function for closed strings 
has been calculated using BPS invariants. The open string story is more complicated. 
Since the A model geometry is not a cotangent bundle the Chern-Simons theory receives 
open string instanton corrections [7]. Nevertheless, the corresponding CS free energy 
has been calculated [H] and the Kahler parameters for open and closed strings have been 
related. An interesting half integral shift in the second Kahler parameter is observed. 
We want to stress that since the shift is not seen at the tree level the only available 
calculation that confirms this shift is that of Diaconescu, Florea and Grassi. A more 
simple alternative calculation would provide us with an additional insight. 

Fortunately, there is a way to write the CS partition function as a matrix model 
integral [HI IH] • In addition, the mirror B model geometry can be obtained by taking 
the large N limit of the matrix model. The mirror to the local conifold is known, we 
will show that the spectral curve of the corresponding matrix model matches this B 
model geometry, thus providing further evidence of the proposed duality. Moreover, 
the dictionary between open and closed string sides is derived using not only the free 
energy but also geometry, importantly we confirm the half integral shift found in [H]. 
This is the main result of the article. We also calculate the genus zero free energy using 
the matrix model. The technique we use to get the spectral curve is inspired by [TU] 
where spectral curves for hermitian two matrix models are studied. 

In the second section we give a concise description of the CY geometry from the A 
and B model sides. The third section is devoted to the matrix model integral and the 
equations of motion. The spectral curve is found and discussed in section 4. 
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2 Geometry 



In this section we briefly introduce the underlying CY geometry. For details see |Slini- 
We start with the deformed local conifold Y in the A model. The toric diagram for 
Y has an extra vertex compare to the toric diagram of the deformed conifold. N 
D-branes are wrapped on the lagrangian S"^ in this geometry. The topological open 
string A model is U{N) CS on with open string instanton corrections [HI that have 
an interpretation of non local Wilson loop effects in the CS theory ITj. The size of the 
is a complex structure moduli therefore A model amplitudes are insensitive to it. 
Taking the to the zero size leads to a singular CY that locally looks like a conifold. 
The singularity can be resolved by blowing up a P^. The geometric transition maps the 
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Figure 1: The web of large N dualities and mirror symmetries. 



open topological model on the deformed local conifold Y to the closed string model on 
the resolved CY Y that has two Kahler parameters, t and s, the sizes of the two P^'s. 
In terms of the gauge theory the transition leads to the large N limit of the CS theory. 
If embedded in the full string theory this geometric transition describes dissolving N 
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D6-branes wrapped on the vanishing and replacing them by a flux. This large N 
duality was conjecture by Aganagic and Vafa in [5^ and further studied by Diaconescu, 
Florea and Grassi in |6'. 

There is the mirror side of the story. The mirror to Y is the resolution of the 
following singular CY |1H IT^ 

a;^ = (e"-* - l)(e" + e^ - 1). (2.1) 

The singularity gets resolved by blowing up a holomorphic at the points where the 
r.h.s. of equation 1)2.11) is zero. Let denote this CY as X. The branes mirror to the A 
model branes wrap this P^. The effective topological open string theory is holomorphic 
Chern-Simons [7j. The large N transition for this geometry closes the web of mirrors 
and large N transitions giving us the mirror to Y, the A model resolved CY. The mirror 
of Y is the deformed CY X 

xz = {l- e-")(l - e"-*) - e-\l - e""*"^) (2.2) 

which is written in the fiat coordinates [TH]. The web of the mirror maps and the large 
N transitions is summarised on fig^ 

The genus zero s dependent part of the free energy for closed strings in the resolved 
geometry Y was found in j3] by computing classical periods in the mirror B model ge- 
ometry X. It was compared with the disk computations for open strings. Furthermore, 
the s dependent part of the all genus partition function for closed strings was derived 
using BPS invariants. The t dependent part of the all genus free energy was obtain in 
jH] by compactifying the model and calculating BPS invariants. The open string com- 
putations are more complicated since they involve the calculation of non local Wilson 
loops in the gauge theory. The open string partition function was obtained in [H] and 
matched to the closed string partition function. The dictionary between the open and 
closed string parameters found by Diaconescu, Florea and Grassi, 

s = -gsN i = top + (2.3) 

contains the curious half integral shift in the second Kahler parameter i. This shift 
is not seen at the tree level. The shift is interpreted as a contribution of degenerate 
open string instantons. So far this is the only available open string computation that 
confirms this shift. 



4 



An alternative computation of the CS partition function can be done by using 
matrix models. Moreover, the spectral curve of the matrix model is a part of the mirror 
to the large N dual geometry. In this context, CS matrix models are a powerful tool for 
studying relationship between closed and open string parameters. The corresponding 
matrix models were introduced in jU] . They differ from more familiar hermitian models 
by the measure factor. The integration goes over the unitary group. CS matrix models 
have been used to study the large N dualities for the conifold 0, the orbifold of the 
conifold [Hj and the lens spaces [13 HE] ■ 

In this paper we find the spectral curve of the CS matrix model for the local conifold 
and confirm eq. (j2.3|) . 



3 Matrix model 

The local conifold under consideration contains the lagrangian over which N D-branes 
are wrapped. Since the geometry is not a cotangent bundle the CS theory gets modified 
by open string instanton corrections [7]. The corresponding two matrix model for the 
CS theory has been proposed in ^ and reads 

Z - JY[ du,dviA{u)A{v)exp (-^ J^'^UiV, + Vi{vi) + V2{ui)) \ , (3.1) 



A(x) = n2sinh(^^ (3.2) 



where the measure factor is 

2sinh(^'~^' 

i<j 

and the potentials are given by 

V,iv)=tv, (3.3) 
^^(«) = -E;^^""- (3-4) 

n=l 

The integration measure is the Haar measure on the space of two commuting unitary 
matrixes. The two extrema of the matrix model action correspond to the north and 
south poles of the holomorphic in the blown up B model geometry ()2.ip . That is a 
characteristic feature of CS matrix models, their potentials encode the mirror B model 
geometry. It has a simple explanation since CS matrix models can also be seen as 
coming from holomorphic CS on the B model side. 
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The matrix model has two sets of eigenvalues. They interact only through the 
linear term ^^UiVi. This gives a hint that one should look at a similar hermitian two 
matrix model which has two matrixes X and Y that interact though TiXY in the 
matrix model action. There is a vast literature on the subject, see for example |17j . 
The one of a particular interest for us is the paper by Kazakov and Marshakov [TU] 
where a simple method for finding spectral curves for hermitian two matrix models is 
proposed. The model also has applications to 2D string theory jTHj . 

Anticipating taking the large N limit denote s = QgN as the 't Hooft parameter. 
Let introduce two resolvents, one for each set of eigenvalues, 

^i(^) = OsY^^ coth I ^ „ ^ ) , (3.5) 



coth ( - 

i ^ 

uJ2{u) = 93^ coth 2 ^ j ' ^^'^^ 

i ^ ' 

They enjoy the following boundary conditions 

lim ci;i^2(a^) = ±s. (3.7) 

The equations of motion (EM) can be conveniently written in terms of the resolvents 
as 

Ui = t+ -iji!i{vi), (3.8) 

=ln(l-e"0 + ^^2(ni). (3.9) 

There is another form of EM that resembles a hermitian two matrix model. Let U = e^, 
V = e^ Wi{V) = e^i(^)/2 and W2(?7) = e'^2{c/)/2 ^^^^^ 

U = e'Wi{V), (3.10) 



V = {1-U)W2{U). (3.11) 

Our goal is to take the large N limit of the above equations and derive the spectral 
curve. 
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4 Spectral Curve 



In this section we consider the large N hmit of the CS matrix model and find the 
corresponding spectral curve. The spectral curve of the matrix model is a part of the 
large N dual B model geometry. The goal is to compare the spectral curve of the 
matrix model to X, the mirror to the conjectured large N dual Y. 

We want to proceed along the lines of jTUI by Kazakov and Marshakov, where the 
spectral curve for the hermitian two matrix model is derived. Although hermitian 
models are quite different from CS matrix models there is an obvious analogy between 
their EM observed in dTTIH) and (ITTTll . 

In the large N limit the eigenvalue distribution becomes a continuous function and 
poles of the resolvents spread into the cuts. Since U and V in the EM in the large N 
limit lie on some cuts the functions Wi(f/) and W2(y) in (j3.1Uj) . (j3.1H) are understood 
as principal values. The EM determine the behavior of the functions U{V) and V{U) 
at infinity and at zero . In addition, as explained in jTU] these two functions in the 
large N limit are inverse of each other 

UiVm = U. (4.1) 

Furthermore, they completely fix the spectral curve of the model. Either one of them is 
actually the spectral curve. The functions U {V) and V{U) have the following behavior 
at infinity 

f/(\/) ~e*+^/2 + 0(l/y), (4.2) 

V{U) ~ + 0(t/°) (4.3) 

and at zero 

U{V) ^ e'-'/^ + 0{V), (4.4) 

V{U) ^ e-'/^ + 0{U). (4.5) 

To find the spectral curve we are going to apply the method proposed by Kazakov and 
Marshakov for hermitian two matrix models. It turns out that our case when U{V) 
goes to a constant at infinity is more subtle then the general one. So we start with the 
general case and briefiy review the technique. We suppose that U{V) ~ for big V 
and V{U) ~ f/™' for big U. Then the spectral curve is given by 

F{U, V) = + - VU"^ + ... = 0, (4.6) 
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where the function F{U, V) is polynomial and only the first three terms with the highest 
powers of U and V are written explicitly. The functional form of F{U, V) is completely 
fixed by the given behavior at infinity. The coefficients of the lower power terms can 
be found using the known behavior of U{V) and V{U) at zero and by fixing certain 
period integrals of the spectral curve. Note that if, for example, n = 0, terms of the 
form [/^ — U^"^"^ + in F{U^ V) for any integers k and / are completely consistent 
with the given asymptotic at infinity. It is clear that for this degenerate case one can 
not determine the functional form of the spectral curve just by looking at the functions 
U{y) and V{U) at infinity. However, it is reasonable to consider the matrix model 
(j3.H) as a smooth limit of a matrix model with a potential that has higher powers of 
V . Let Vi(f) be such that 

Vi{v) = In (e* - e^V) . (4.7) 

When r — oo the above potential reduces to Vi{v). With this new potential 
U{y) ~ V"^ at infinity and so the formula (j4.6|) with n = 2 and m = 1 can be ap- 
phed leading to 

F{U, V) = aV^ + (3U^ + ^V'^U + SV"^ + eUV + eU + CV + 1. (4.8) 

The next step is to see what coefficients survive taking r — oo limit. As it is shown 
in the appendix the spectral curve for the cubic potential looses three terms in this 
limit and takes the form 

F{U, V) = (3U^ + eUV + eU + (¥ + 1 =0. (4.9) 

There are four coefficients that can be found by the known behavior at infinity and at 
zero. Carrying on this calculation one gets 

F{U, V) = U^ + e-'I^UV - e'U - (1 + e'~'/^)V + e*-^/^ _ q_ ^4^0) 

The simplest check of the above expression is to take the limit s = 0. In this limit all 
eigenvalues must be located at the classical extremum of the matrix model action. That 
is indeed true. The spectral curve in the classical limit splits into the two independent 
equations V = and V = 1 — U that are precisely the location of the holomorphic 
over which the branes are wrapped. The spectral curve also gives us the resolvents 

Wi(V) = ^ (^1 + e*-^/2 - e-^/V + ^J (e-'V - 1 - e'-'/'^f + Ae'V - 4e*-^/2^ , 

(4.11) 
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W2(f/) 



(4.12) 



U - e*+V2 ' 

where the plus sign in front of the square root is chosen in order to have the correct 
hmit at infinity. Note that the point v = — oo where the resolvent uji{v) is equal to 
—s is located on the other sheet of the spectral curve which corresponds to the minus 
sign before the square root in ()4.11|) . In the next section we calculate the genus zero 
free energy as an integral over u of the spectral curve so we also need the expression 
for the spectral curve in the u variable 



v{u) = In ■ 



/2) 



(4.13) 



We now compare our spectral curve to the spectral curve coming from CS on or 
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Figure 2: The spectral curve for the CS matrix model. The two contours are the 
noncompact cycles used in the calculation of the genus zero free energy. 

the A model conifold Since the toric diagram for the resolved A model geometry 
has an extra vertex we expect an extra cylindrical sheet in the mirror B model. Let 
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look at the function u{v) = t + \nWi{v). The corresponding Riemann surface consists 
of the three cyhndrical sheets. The physical sheet is the one where the resolvent uji{v) 
is finite at infinity. The physical sheet is glued to the second cylinder along the square 
root cut. On this second sheet the resolvent uji{v) has the minus sign in front of the 
square root and is equal to — s at t> = — oo. The second sheet has the logarithmic cut 
that starts at f = — s/2 and goes to infinity. The logarithmic cut connects the second 
sheet to the third one that does not have any square root cuts. The logarithmic cut is 
responsible for the extra vertex in the A model toric diagram that differs it from the 
conifold. The square root cut shrinks to the zero size when s = while the logarithmic 
cut is present even in the classical limit although the location of the branch point gets 
modified by quantum corrections. The spectral curve is depicted on figI21 

It is easy to see that the spectral curve looks very similar to ()2.2j) that is the mirror 
to the conjectured large N dual of Y. It turns out that there is a simple coordinate 
transformation 

v = -v + s/2 (4.14) 

u= -u (4.15) 
that brings the spectral curve to the form 

^t-s/2^2u+v _ gt+./2g2« e*""/=^)e"^ + + e" = (4.16) 

which is exactly the non trivial part of the deformed B model CY X in (j2.2j) . Now it is 
easy to read off the relationship between the open and closed string Kahler parameters 

s = -s i=-t + s/2, (4.17) 

where i and s are the closed string flat coordinates. Since the open string Kahler 
parameter t that is used in the matrix model is equal to —top of 6J the above equations 
are precisely those found by Diaconescu, Florea and Grassi ()2.3|) . It confirms the large 
N conjecture and the shift by s/2 in the second Kahler parameter. 

Therefore we have derived the mirror to the large N dual of the local conifold. It 
coincides with one proposed by Aganagic and Vafa in (5]. This is a strong evidence 
that the large N conjecture is indeed true. 

To further study the large N transition one has to find the matrix model free energy 
and compare it with the closed string calculation. 
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5 Free energy 

In this section we calculate the genus zero free energy. As for hermitian matrix models 
the derivatives of the free energy with respect to t and s can be written as integrals of 
the resolvents over some noncompact cycles. A useful relationship that we are going 
to explore through this section is [T3] 

^ In sinh ^ 2^' ) ^ uj{z)dz - ^^Xi, (5.1) 

where A is a point at infinity on the physical sheet of the resolvent and x is a point 
on the cut. All A dependent terms are omitted. The sum g.sYli'^i is very important 
and represents a period integral, meanwhile the sum Qs Y2i '^i — is trivial and will 
be omitted. Since we have to sum over both u and v eigenvalues it is convenient to 
introduce A„ and A^, points on the physical sheet at u and v infinities respectively. 

It immediately follows that the derivative of the genus zero free energy Fq with 
respect to t can be cast into the form of the integral of the spectral curve 

dtFo = gsY^Vi = gsJ^H^ = - (5.2) 

i i 

Here, the summed over i EM have been used and all polynomial in t and s terms are 
skipped. v{u) is given by (j4.13|) . The integral can be easily evaluated leading to 

oo ^ 

d,Yo = Y.—2 (e"^*-^/'^ - e<'+''^^) . (5.3) 

n=l ^ 

The derivative of the free energy with respect to the 't Hooft parameter is the variation 
of the matrix model action due to eigenvalues u and v{u) added to the square root cut 



dsFo = -uv + Vi{v) + V2{u) + J] In U sinh j + (7s J] In U sinh j . 

(5.4) 

The two last sums in the above equation can be converted to the integrals of the 
resolvents by means of ()5.H1 . Using ()5.2j) the expression can be further transformed to 
the form 

Ai, /"All 1 /"Au 



SsFq = —UV — I u{v)dv — I v{u)du ^ I v{u)du. (5.5) 
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The last step is to use that u and v are related by the spectral curve (|4.1Up . This takes 
us to (see figI21) 

dgFo = — / v{u)du + - / v{u)du (5.6) 

Js/2+t 2 Jq 

with the result of the integration being 

dsFo = -f2-2 + ^e"(*-^/2) + le"(*+^/2)) . (5.7) 

n=l ^ V / 

The energy becomes 

n=l 

With the identifications ()4.17j) this free energy matches the closed string free energy 
found in 0. 

The fact that dgFo is the combination of the integrals over the two independent cy- 
cles is a consequence of the half integral shift in ()4.17|) . The open string flat coordinate 
t is not a flat coordinate for the closed string model. In the matrix model language 
if we interpret the matrix eigenvalues as D-branes the shift occurs because the anti 
D-branes at m = form a bound state with the D-branes spread between u = t — s/2 
and u = t + s/2. Moving two D-branes to inflnity must be supplemented with pulling 
one anti D-brane to inflnity. 

It leads to another check of the large N transition. The free energy calculated from 
the open string prospectives is the same as the closed string energy. 

Conclusion 

We have provided a strong evidence that the conjectured large N transition is indeed 
true. We see that the second Kahler parameter t does receive the half integral shift 
that supposedly comes from the degenerate open string instantons. 

We have found the method of calculating spectral curves for the wide class of 
two matrix Chern-Simons models. It is easy to see that the method works for any 
polynomial (after exponentiating ) potentials. The spectral curve is given by (j4.6p . 

It is interesting to consider more general Chern-Simons matrix models, for example, 
the matrix models for the geometries considered in [T^. Given a relative simplicity 
of matrix model calculations it is conceivable that Chern-Simons matrix models can 
provide many new interesting examples of large N transitions. 
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Another intriguing question is to relate universal correlators for 3D partitions in 
CY crystals to CS matrix models pUj . 
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A From cubic to linear potential 

The new EM with the cubic potential are 

U{V) = (e* - e'V^)Wi{V), (A.l) 

ViU) = {l-U)W2{U), (A.2) 

where parameter r is intended to be taken to minus infinity at the end. Our concern 
is only the behavior of U{V) and V{U) at infinity which is given by 

U{V) ~ -e'^e^/V^ - e^'AV + {e'+'/^ - e'B) + 0{1/V), (A.3) 

\/(f/) ~ -t/e'^/2^0(f/°). (A.4) 

Here A and B stand for the expansion coefficients of the resolvent Wi {V) . The resolvent 
has to be expanded up to l/V"^ in order to keep the constant term in ()A.3|) . The later 
is needed since after putting to zero eq. ()A.3|) must reduce to ()4.2|) . The spectral 
curve that is consistent with the above behavior at infinity is of the form ^U] 

aV^ + (3U^ + -fV^U + 6V^ + eUV + eU + + 1 = 0. (A.5) 

We want to find what coefficients in the above equations are zero if is zero. As V 
goes to infinity we have to zero all terms of order V^, and V'^ which gives 

e'+'/^(3 -7 = 0, (A.6) 
a + 2pe^'^'/^A - -fe'A - ee'^+^Z^ _ g, (A.7) 

e^'A^p - 2/5(e*+'/2 - e"S)e"+'/2 ^ ^(e*+^/2 _ g'-^) + 6- ee'A = 0. (A.8) 
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As U is taken to infinity only terms of order liave to be kept 

-ae"/2 + 7 = 0. (A.9) 

Clearly, when = one must have a = 7 = 5 = 0. Therefore, the appropriate 
functional form of the spectral curve with the linear potential is 

(3U^ + eUV + eU + CV + 1 = 0. (A.IO) 



References 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 



[9 



R. Gopakumar and C. Vafa, "M-theory and topological strings. I," 
iarXiv:hep-th/9809187| 

R. Gopakumar and C. Vafa, "On the gauge theory/geometry correspondence," 
Adv. Theor. Math. Phys. 3, 1415 (1999) |arXiv:hep-th/9811131 . 



H. Ooguri and C. Vafa, "Worldsheet derivation of a large N duality," Nucl. Phys. 
B 641, 3 (2002) |arXiv:hep-th/ 0205297|. 

T. Okuda and H. Ooguri, "D Branes and Phases on String Worldsheet", 
|arXiv:hep-th/0404101| 

M. Aganagic and C. Vafa, "6*2 Manifolds, Mirror Symmetry and Geometric En- 
gineering", arXiv:hep-th/0110171 



D. -E. Diaconescu, B. Florea, A. Grassi, "Geometric Transitions and Open String 
Instantons", Adv. Theor. Math. Phys. 6 619 (2003), arXiv:hep-th/0205234| 

E. Witten, "Chern-Simons gauge theory as a string theory," Prog. Math. 133, 637 



(1995) |arXiv:hep-th/9207094 . 



M. Marino, "Chern-Simons theory, matrix integrals, and perturbative three- 
manifold invariants," |arXiv:hep-th/0207096, 

M. Aganagic, A. Klemm, M. Marino and C. Vafa, "Matrix model as a mirror of 



Chern-Simons theory", JHEP 0402 010 (2004), |arXiv:hep-th/0211098 . 



14 



[10] V. Kazakov and A. Marshakov, "Complex Curve of the Two Matrix Model and 



its Tau-function" , J.Phys. A36 3107 (2003), arXiv:hep-tli/0211236 



K. Hori and C. Vafa, "Mirror symmetry," |arXiv:hep-tli/0002222 



K. Hori, A. Iqbal and C. Vafa, "D-Branes And Mirror Symmetry", 
|arXiv:hep-tli/0005247| 

M. Aganagic, A. Klemm and C. Vafa, "Disk Instantons, Mirror Symmetry and 
the Duality Web", Z.Naturforsch. A57 1 (2002) , | arX iv:hep-th/0105045| 

N. Halmagyi and V. Yasnov, "Chern-Simons matrix models and unoriented 
strings" , arXiv:hep-th/0305134'| 

N. Halmagyi and V. Yasnov, "The spectral curve of the lens space matrix model" , 



|arXiv :hep-t h/0311117| 

N. Halmagyi, T. Okuda and V. Yasnov, "Large N Duality, Lens Spaces and the 
Chern-Simons Matrix Model", JHEP 0404 014 (2004) , [a^v:hep-th/0312l45 . 



[11 
[12 

[13 

[14 

[15 

[16 

[17] J. Alfaro, "The Large N Limit of the Two Hermitian Matrix Model by the Hidden 
BRST Method", Phys.Rev. D47 4714 (1993), |arXiv:hep-th 7920709T ; M. Stau- 
dacher, "Combinatorial Solution of the Two-Matrix Model", Phys.Lett. B305 332 
(1993),|arXiv:hep-th/9301038|; P. Zinn- Justin, "Universality of correlation func- 
tions of hermitian random matrices in an external field", Commun. Math. Phys. 
194, 631 (1998), |arXiv:cond-mat /97050441 ; B. Eynard, "Large N expansion of 
the 2-matrix model", JHEP 0301 051 (2003), |arXiv:hep-th/ 0210047| 

[18] S. Alexandrov, V. Kazakov and 1. Kostov, "2D String Theory as Normal Ma- 
trix Model", Nucl. Phys. B667 90 (2003), arXiv:hep-th/0302106.; S. Alexan- 
drov, V. Kazakov and D. Kutasov "Non-Perturbative Effects in Matrix Mod- 
els and D-branes", JHEP 0309 057 (2003), |arXiv:he p-th/0 306177| ; V. Kazakov 
and I. Kostov, "Instantons in Non-Critical strings from the Two-Matrix Model" , 
arXiv:hep-th/0403152 

[19] D.-E. Diaconescu, B. Florea and A. Grassi, "Geometric Transitions, del Pezzo 
Surfaces and Open String Instantons", Adv. Theor. Math. Phys. 6 643 (2003), 
|arXiv:hep-th/0206163] 



15 



[20] R. Dijkgraaf, A. Sinkovics, M. Temurhan "Universal Correlators from Geometry" , 
|arXiv:hep^tli/0406247| 



16 



